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ABSTRACT 

We present a method to estimate and map the two-dimensional distribution of dust extinction in 
the late-type spiral galaxy NGC 959 from the theoretical and observed flux ratio of optical V and 
mid-IR (MIR) 3.6 /im images. Our method is applicable to both young and old stellar populations 
for a range of metallicities, and is not restricted to lines-of-sight toward star-formation (SF) regions. 
We explore this method using a pixel-based analysis on images of NGC 959 obtained in the U-band 
at the Vatican Advanced Technology Telescope (VATT) and at 3.6 /xm (L-band) with Spitzer /IRAC 
We present the original and extinction corrected GALEX far-UV (FUV) and near-UV (NUV) images, 
as well as optical UBVR images of NGC 959. While the dust lanes are not clearly evident at GALEX 
resolution, our dust map clearly traces the dust that can be seen silhouetted against the galaxy's 
disk in the high-resolution HST images of NGC 959. The advantages of our method are: (1) it 
only depends on two relatively common broadband images in the optical U-band and in the MIR 
at 3.6 /im (but adding a near-UV band improves its fidelity); and (2) it is able to map the two- 
dimensional spatial distribution of dust within a galaxy. This powerful tool could be used to measure 
the detailed distribution of dust extinction within higher redshift galaxies to be observed with, e.g., 
the fiS'T/WFC3 (optical-near-IR) and JWST (mid-IR), and to distinguish properties of dust within 
galaxy bulges, spiral arms, and inter-arm regions. 

Subject headings: dust, extinction — galaxies: individual (NGC 959) — galaxies: spiral — galaxies: 
structure 



1. INTRODUCTION 

Dust extinction is a longstanding problem when study- 
ing stellar populations within our Galaxy and in ex- 
tragalactic obj ects (e .g., iTrumpleH 11930': 'Mathis et al.l 
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Roussel et al. 



1982; Caplan & Deharveng I198E 
river et al.l 120081) . Variations 
the spatial distribution of in- 
different effects on the light 



in the amount and 

terstellar dust have 

from background stell ar pop ulations ( e.g.. lElmegreenI 



1991 



1980; Wahcrbos & K ennicuttI 1^ )88; W aller et al.l — . . 
Witt et all 119921 : iCalzetti et al.lTl994t iDeo et al.l l2006n . 
Measuring and correcting for dust extinction in indi- 
vidual galaxies is important to understand the true 
nature of their stellar populations, especially when 
one aims to study small-scale structures such as Hii 
regions of spiral galaxies, which ten d to be par- 
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ticularly dusty (e . g., [P cter 
iScoville et aPlloOlUCaSetti et al.ll200 5). 

Many different methods are used to estimate 
dust extinction within target galaxies. Some com- 
monly used methods involve: (1) the ratios of 
Hydrogen recombination- line fluxes for Hii re gions . 
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UV fluxes (e.g., IBuat fc Xul Il996l : ICalzetti et aD l2000L 
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lej^ lCalzetti et all 11991 iMeurer et all 119991: [Selil 
l2002bf) : and (4) the CO column density from (sub- 
)millimeter observations of, e.g., ^^CO(J=1-0) line 
and J = 1-0 and J = 2 - 1 transitions of^ ^ ^CO and C^^O 
(lEncrenaz et all 119751: i Dickman' '1 9781: iFrerking et alJ 
198a iHa riunpaa fc Mattila.l99& : .Havakawa et al.lll999L 
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Even though all these methods measure dust ex- 
tinction in some way, each has its limitations. The 
Balmer decrement, Ha/H/3, and the Ha/Paa flux ra- 
tios are generally limited to lines-of-sight toward Hii 
regions. These Hii regions can be distributed all 
over a galaxy, but cover only a small fract ion of an 
entire galaxy disk (e.g., IScoville et al.l 120011 ). Accu- 
rate measurements of the underlying Balmer absorp- 
tions from moderate-resolut ion sp ectra are also r equire d 
(|Petersen fc GammelgaaTdl 119971 : IBoissier et all 12004"). 
The FIR/UV flux ratio can be used throughout a galaxy, 
although the spatial resolution in the FIR is gener- 
ally poor: about ~ 20" for t he Mid course Space Ex- 
periment (MSX) at 4.2-26 /xm (iPrice et al. 200|i^ ~ 40" 
for Spitzer /MIPS at 160 ^m (Ricke et al. 2003). and 
even l arger (~ 100") for IR AS at 100 ^m f Xu fc HelotJ 
Il996t IBoissier et al.l 120041 ). This becomes a signifl- 
cant limitation in studying small-scale structural fea- 
tures, and limits access to only the very nearest galax- 
ies. Unlike the FIR images, the UV images can have 
much higher resolut ion (FWHM< 0"1) with the HST 
WFP C2 and ACS (jTrauger et all ir994: Sirianni et al 



2005), or FWHM<4-6" for GALEX (Morrissev et al 



20051) . A major complication of methods involving UV 
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filters (FIR/UV flux ratio and UV spectral slope) is 
that the UV emission is significantly affected not only 
by dust, but also by the age and metallicity of stel- 
lar populations. Finally, sub-millimeter and millime- 
ter observations are used to measure the column den- 
sity of CO molecules directly, and hence — under cer- 
tain assumptions — that of molecular hydrogen, H2. 
The CO (and hydrogen) column density is converted 
to extinction using the correlation found bet ween CO 
and dust-extin ction in nearby galaxies (e.g., 'Dickman' 
Il978; Bachiller fc Cerni charo 1986; Komugi ct al. 2005). 
The typical resolution at mi llimeter wavelengths is 
poor (FWHM-20"- 2f6; Encre naz et"allll975t iDickmanI 
19781 iFrerking et all [l982: Ha riunpaa k Mattilal Il99a 



Havakawa et al.lll999U2 001: Ha rjunpaa et all l2004f) . "al- 
though ALMA ([Brown et al. 200^) will so on greatly im- 



prove on this. Studies such as that of iBoissier et al.l 
(|2007f l also show that the relationship between the dust 
column density and the amount of extinction is compli- 
cated, and that further study is necessary to derive the 
relationship between the two valu es. 

Without using these methods, IReganI (2000^ investi- 
gated dust extinction using the more commonly used 
ground-based optical BVRI and near-IR (NIR) JHK fil- 
ters to obtain various color maps. The advantage of using 
optical-NIR filters is that the extinction measuremen t 
is independent of the dust temperature (|Reganl I2000D , 
and that the images can have much higher spatial res- 
olutio n at these wavelengths (e.g., Reg an k, Mulcha"evl 
I19991) . Combining theoretical models and color maps 
from optical-NIR images, the spatia l dust distribution 
can then be derived. Fig. 7 of IReganI (|2000l ) shows that 
this method can reconstruct e xtinction-free galaxy im- 
ages. However, as lReganl (|2000l ) points out, there are sev- 
eral issues with this particular method. The first is that 
all of the bandpasses involved are affected by the dust to 
some degree: no single optical-NIR filter directly maps 
either stellar or dust morphology by itself. The second 
issue is that, while the intrinsic colors of the underlying 
stellar populations are well known for the older stellar 
populations, the colors of younger stellar populations in 
star-formation (SF) regions vary, dependi ng on their en - 
vironment and specific properties. Since IReganI (I2OOOI ) 
treated all stellar populations as old stellar populations, 
the extinction measurements for stellar populations in 
SF-regions are less accurate and should be treated with 
care, as their paper points out. 

In this paper, we will map the two-dimensional spa- 
tial distribution of dust extinction in NGC959, using 
images from commonly used optical-MIR filters. Many 
images observed with optical-MIR filters from both 
ground and space are now readily available for large 
samples of galaxies though publicly available archives. 
We will also treat the younger stellar populations sep- 
arately from older stellar populations to measure dust 
extinction in SF-regions. We initially combine data 
spanning the GALEX FUV and NUV, the ground- 
based optical UBVR fihers, the 2MASS NIR JHKs fil- 
ters, and the Spitzer /lYiKC 3.6, 4.5, 5.8, and 8.0/^m 
(MIR) filters. The IRAC 3.6 and 4.5 //m filters are 
not nearly as affected by dust extinction as the NIR 
filters, and therefore are commonly used to trace the 
distribution of stellar populations (e.g., iWillner et al.l 
l2004i) . We use the library of spectral energy distribu- 



tions (SEDs) bv I Anders fc Fritze-von AlvenslebenI (|2003l . 
hereafter AF03), to study possible combinations of filters 
and colors in the analysi s of dust ex tinction. U sing a 
pixel - based analysis (e.g.. IBothun|[T986 ; Abraha m et al.l 
11999; Eskridge et alj 120031 : ILanvon-Foster et all 120071 : 
Wclikala ct al. 2008^ on NGC959, we will demonstrate 
that our method can reveal a two-dimensional distribu- 
tion of dust extinction. 

This paper is organized as follows. In §2, we present 
our method to estimate dust extinction using the SED 
models. In §3 and §4, we test and apply this method to 
all available images for the nearby late-type spiral galaxy 
NGC959. We give a discussion of our results in §5, and 
present our conclusions in §6. 

2. MODELS AND CONCEPT 

Before we analyze any observed images, we first evalu- 
ate theoretical SED models for single stellar populations 
(SSPs) at different ages and metallicities. In this section, 
we describe the SED library that we used in our study, 
and then select the optimal filters for the subsequent SED 
analysis. Once the filters are selected, we describe how to 
estimate the dust extinction, using the theoretical SEDs 
and observed images, through fiux ratios of an optimal 
set of filters. 

2.1. Simple Stellar Population Models 

Among published SSP SED librarie s (e.g., AF03; 
IBruzual fc Charlotl[2003t IMarastonll2005h . we elected to 
use the SED library by AF03. While other SED li- 
braries do not contain information about emission lines, 
this library includes both spectral and gaseous emission 
for young stellar populations. AFG3 has created mul- 
tiple sets of SED lib raries^ of SSPs , using the Padova 
(iBertelli et al.l I1994D and Gen eva (iLeieune fc Schaerei 

misochrone models with IScalol (119861) . ISaloeted 
, and iKroupal (|2002|) Initial Mass Functions 
(IMFs). These SED libraries contain models for metal- 
licities of Z = 0.0004, 0.004, 0.008, 0.02 (Solar), and 0.04. 
Both isochrone models include the thermally pulsing 
asymptotic giant branch (TP-AGB) phase and model 
atmosphere spectra from ILeieune et al.l (|1997[ ). The 
main difference between the two libraries with different 
isochrone models is the time resolution. At. While the 
library with the Padova isochrones have an age coverage 
from 4 Myr to 14 Gyr with = 4 Myr up to an age of 
2.35 Gyr, and At = 20 Myr for older ages, the library with 
the Geneva isochrones have variable time steps. Among 
the available models, we will focus on the ones computed 
with Padova isochrones and the Scalo IMF — adopting 
an upper mass limit of ^ 50 Mq for super-solar metal- 
licity, and ~ 70 M© for all other metallicities. We refer 
the reader to AF03 and references therein for a detailed 
description of this library. The SEDs of SSPs change 
rapidly at all wavelengths, but especially so in the UV 
for ages younger than 100 Myr (see Fig. [T|). Once a 
stellar population reaches the age of ~ 500 Myr, the rate 
of change in SED diminishes with increasing age. As a 



Since the publication of lAnders &: Fritze-von Alvensleberil 
II2003I) . there have been several updates and additions to their 
models, resulting in slightly different combinations of IMFs and 
metallicities from the original description. We used the version 
dated 28 October, 2007. 
http : //www. galev . org/ 
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Fig. 1. — Top panel: T heoretical SEDs from the library of 
lAnders &: Fritze-von Alvens lcbcn (2003). Among the different sets 
of isochrones and IMFs, SEDs for Padova isochrones assume a Scalo 
IMF with metallicity of 0.0004 (solid curves) and 0.04 (dash-dot 
curves) at 6 different ages as indicated. The age-metallicity degen- 
eracy is clearly visible for wavelengths shortward of 4000 A. On the 
other hand, SEDs at A > 6000 A are less affected by metallicity, ex- 
cept for extremely young (t < 100 Myr) stellar populations. Bottom 
panel: Total throughput curves, T(A), for different telescope-filter 
combinations. Throughput for the GALEX filters are scaled up 
by a factor of 10 for better visibility. Vertical dotted lines indicate 
the wavelength coverage for each telescope-filter combination. 



practical subset for our analysis, we consider 10 differ- 
ent ages (t = 4, 8, 12, 52, 100, 500 Myr, 1, 5, 10, and 
13.5 Gyr) for each metallicity, resulting in a model grid 
of 50 different SEDs. In the main panel of Fig. [1] we plot 
the SEDs for six ages and two metallicities that span the 
full range of our adopted model grid. 

2.2. Optimal Filter Selection 

The bottom panel of Fig. [1] shows the total through- 
put curves, T(A), for various FUV-MIR filters consid- 
ered in this study, as well as in the subsequent multi- 
wavelength studies on NGC 959 and on a sample of 45 
galaxies (Tamura et al. 2009b, c, in preparation). Total 
th roughput curves for t he GALEX filters — published 
in lMorrissev et al.] (|2005i ) and on the GALEX web page^ 
— are multiplied by a factor of 10 for better visibility. 
The throughput curves for the VATT filters are the to- 
tal responses of the filters, CCD, and telescope^. The 
throughput curves for the 2MASS JHKs filters are the 
total responses of the filters, CCD, telescope, and atmo- 

* Goddard Space Flight Center, the GALEX Post Launch 
Response Curve Data: 

http : //galexgi . gsf c .nasa.gov/docs/galex/Documents/ 
PostLaunchResponseCurveData .html 

^ Instrumentation for VATT: 
http : //vaticanobservatory . org/vattinst .html 
The response curve used in this paper is for the VATT 2k CCD. 



sphere^. Finally, for the Spiteer /IRAC filters, the total 
throughput curves are calculated using all components 
of the Spitzer instrument''. The vertical dotted lines in 
Fig.[T]indicate the wavelength ranges covered by different 
telescopes and instruments. 

The largest change in SEDs with the increase in age 
is the reduction of flux for wavelengths shortward of 
4000 A. Especially the UV flux — as measured by, e.g., 
the GALEX FUV and NUV filters — decreases by up to 
~ 6 dex relative to the maximum flux level, as an SSP 
ages from 4 Myr to 13.5 Gyr. SEDs are also significantly 
affected by metallicity. For SSPs older than 100 Myr, 
the metallicity mainly affects the wavelengths < 4000 A, 
where the fiux decreases by up to ~ 2 dex from Z ~ 0.0004 
to Z = 0.04. Since our primary goal is to measure dust 
extinction and to map its spatial distribution, filters cov- 
ering wavelengths shorter than 4000 A should be avoided 
to minimize effects from the age-dust-metallicity degen- 
eracy. 

According to AF03, the model SEDs become uncer- 
tain longward of 5 /im. Also, the observed stellar 
continuum at MIR wavelengths redward of 5 /xm is 
contaminated by e mission from polycycli c aromatic hy- 
drocarbons (PAHs: iLeger fc Pugetl Il984f ) and silicates 
(jWillner et al.ll2004D . We thus exclude Spitzer /IRAC 
4.5, 5.8, and 8.0 /im filters. On the other hand, the effects 
from dust extinction and emission by PAHs and silicates 
reach a minimum near 3 .5 /xm (L-band; e.g., F azio et al.l 
l2Oi0llWillner et al.ll200l . The IRAC 3.6 ^m filter there- 
fore provi des the most reliable stellar population tracer 
(see, e.g.. lKennicutt et al.|[200l iHelou et al.ll2004D . 

From the filters with 0.4<A<3.5^m — optical VR 
and NIR JHKs in our study — we need to select (at 
least) one more filter to trace dust extinction. Since dust 
extinction is much sma ller in the NIR than in the optical 
V- and j?-bands (e.g.. Ivan Houtenlll961t ICardelh et al.l 
[1981 ICalzetti et al.lfigM IGordon et al.ll2003D . NIR fil- 
ters are not optimal for this purpose. Another reason 
to avoid ground-based NIR observations is absorption 
due to atmospheric water vapor. The amount and un- 
certainty due to this absorption depend on atmospheric 
condi tions, as well as the locati on of the observations 
(e.g.. iNitscheh^ 119881: ICohen et al. 2003). Dust extinc- 
tion is stro nger in the V-band (e.g., Calz etti et al.lfl99^ 
IGordon eTa l. 2003), while the metallicity effects are 
slightly weaker (see Fig. 1). We therefore choose the V 
and 3.6 /im filters to globally trace the dust extinction. 
In the following section, we will explore the theoretical 
y-to-3.6 /im fiux ratio in detail, and explain how we can 
use this ratio to estimate dust extinction. For other stud- 
ies, in general, it should be noted that filters like SDSS 
5, HST F555W, F550M, and F606W would be adequate 
substitutes for y-band after calibration of the theoretical 
flux ratios. 



2.3. Theoretical V-to-3.6 Flux Ratio Map 

^ 2MASS All-Sky Data Release Explanatory Supplement 
Facilities and Operations: 

http : //www . ipac . caltech. edu/2mass/releases/allsky/doc/ 
sec3_lbl .html 

^ Spitzer Science Center: IRAC: Spectral Response: 
http : //ssc . spitzer . caltech. edu/ irac/spectral_response .html 
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Fig. 2. — Map of the theoretical y-to-3.6^ m flux ratio, 0v, 
calcul ated from an SSP library of Anders & Fritz e-von AlvenslebenI 
II2003I) for varying age and metallicity. Contours are plotted at 
f3v =0.5 (solid curves) and 2.0 (dashed curve). The orange stippled 
region (and the orange band in the color bar) indicates the regions 
in the metallicity-age space where the flux ratio is 0.75 < f^y < 1.50. 
While f^Y varies rapidly for young — t < 500 Myr or log (tMyr)jS2.7 
— stellar populations, it is relatively stable and insensitive to age 
and metallicity for older stellar populations. 

Given the SEDs from the AF03 hbrary and the 
throughput curves of the optical V and IRAC 3.6 /ini 
filters (see Fig. [ij, we explore the range in the theo- 
retical flux ratio, (3v = fv / fs.e fj.m, as a function of age 
and metallicity. The resulting flux ratio map in age- 
metallicity space is shown in Fig. [2l From the input of 
50 SEDs described above, both age and metallicity are 
resampled logarithmically using 100 steps with Spline in- 
terpolation to map the variation in f3v values over the full 
parameter ranges. The minimum ratio, /3y_min = 0.39, 
occurs for the most metal-rich (Z = 0.04 = 2Zq) and 
oldest {t = 13.5 Gyr) SSP, while the maximum ratio, 
/Jy.max = 17.4, is found for the most metal-rich and the 
youngest {t = A Myr) SSP. Gontours for /3v = 0.5 (white 
solid curves) and 2.0 (white dashed curve) are drawn 
in Fig. O A band of orange dots — also included as 
an orange band in the color bar — indicates the re- 
gion where the V-to-3.6 fim flux ratio spans the range 
0.75 </3v< 1.50. 

For stellar populations older than 500 Myr, or 
log(tMyr) ^ 2.70, the theoretical flux ratio changes little 
as a function of age and metallicity. Fig. [2 shows that 
Pv ranges mostly between ~ 0.5 and ^ 2.0 . On the other 
hand, for stellar populations younger than 500 Myr, /3y 
changes more dramatically, depending on their age and 
metallicity. Especially for SSPs with super-solar metal- 
licity, (3v changes from ~ 17 (red) for the youngest SSPs 
to <0.5 (dark purple/black) in a matter of only a few 
100 Myr. For the young SSPs with sub-solar metallici- 
ties, the value of /3y is more stable, but still changes much 
more rapidly than for SSPs older than 500 Myr. Among 



these young sub-solar metallicity SSPs, Fig.[2]shows two 
distinctive ranges of (3v values: (1) a band color-coded in 
green (9 < /3y < 13) for SSPs younger than 10-100 Myr; 
and (2) a band coded in blue (4 < /3y < 7) for SSPs up 
to - 500 Myr. 

The distribution of 0v values in Fig. [2] confirms the 
assumption of, e.g.. iReganl (|200Clf ): that old stellar pop- 
ulations have relatively constant color, while young stel- 
lar populations change their colors more rapidly depend- 
ing on their properties and environments. An important 
corollary of the former is that a mixture of SSPs with 
ages larger than 500 Myr will have a Pv value that is 
very similar to that of a single SSP at those ages. Our 
analysis of the theoretical y-to-3.6 /Ltm flux ratios further- 
more shows that — at least for young stellar populations 
with sub-solar metallicities — the theoretical flux ratio 
can be approximated as being constant for stellar popu- 
lations in either of the two age ranges {t < 10-100 Myr or 
10-100 < i < 500 Myr) described above. In the following 
section, we will derive a method for using these theoreti- 
cal and observed T^-to-3.6 /.tm flux ratios to estimate the 
amount and spatial distribution of dust extinction in the 
y-band. 

2.4. Estimating the V-hand Dust Extinction, Ay 

In the previous section, we showed that the theoretical 
y-to-3.6/im flux ratio, /3y, is well-behaved, particularly 
for sub-solar metallicities. To estimate the dust extinc- 
tion for an observed stellar population, we first have to 
determine its approximate age and metallicity. The ap- 
pearance of horizontal bands of different Pv values in 
Fig. [2] indicates that an age determination is more im- 
portant than the determination of metallicity. Based on 
the Pv values and corresponding age range for different 
metallicities as described above, the accuracy required 
in determining age is ~ 0.5 dex for extremely metal poor 
stellar populations, and 1 dex for sub-solar to solar 
metallicity and older {t > 500 Myr) stellar populations. 
For stellar populations with super-solar metallicities, the 
accuracy required is also ~ 0.5-1 dex for most ages, un- 
less they are extremely young (t<10 Myr). The use of 
color-magnitude diagrams (CMDs) and color-color dia- 
grams, therefore, should suffice to estimate the age of a 
stellar population and select a theoretical Pv value from 
Fig. [21 By comparing the theoretical and observed fiux 
ratios, we can infer the amount of flux missing in the 
V-band observation due to dust. The missing flux is a 
very robust property of the dust distribution, indepen- 
dent of the geometry and flne structure of obscuration. 
The amound of extinction by dust along the line-of-sight 
that this stellar pouplation suffers can then be estimated 
from this missing V-band flux. 

The dust extinction, Ay, is defined as 



Ay = {mv — mvfi) ■. 



(1) 



where my and my.o denote the observed and intrinsic 
(extinction-free) V-band magnitudes, respectively. Even 
though the extinction-free magnitude is not known for an 
observed stellar population, we can use the selected f3v 
value and observed 3.6 /xm flux to estimate the predicted 
mv.o- Since we assume that the obs erved 3.6 jum flux 
is unaffected by dust extinction (e.g.. iFazio et al.l 120041 : 
iWillner et al... 20041 . we can estimate the extinction- free 
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y-band flux by multiplying the selected theoretical (3v 
value with the observed 3.6 /im flux: 

fvfl = PV ^ f 3.6 flTTL ■ (2) 

Therefore, equation ^ can be rewritten as 

Av=mv- [-2.51og(^y x /a.eMm) - V^p] , (3) 

where T4p is the zero-point magnitude for the y-band. 
Compared t o popul ar meth ods such as th e UV spectral 
slopes (e.g.,|BeU||2002b; Ko ng et al.ll2004D and the FIR - 
to-UV flux ratio (e.g., Dale et al.ll2001l : lBuat et al.|[200l . 
this method is much simpler. 

Even though our prescription is simple, there remain 
many possible sources of error. In addition to the usual 
uncertainties due to observational measurements, there 
are several sources of uncertainty in the selected theo- 
retical Pv value. Only when the age and metallicity are 
known for a resolved stellar population can we deter- 
mine /3y for this stellar population with minimal error. 
For unresolved stellar populations, determining the age 
and metallicity becomes much harder as a result of light- 
blending from intermixed and superposed stellar popu- 
lations. As we have shown above, the metallicity depen- 
dency is not as strong as the age dependence in deter- 
mining the likely Pv value. We therefore are still able to 
estimate the Pv value, if an approximate age of the stel- 
lar populations can be determined. Even though stellar 
population composites span a range in age, their luminos- 
ity tends to be dominated by their younger components 
at most wavelengths (see Fig. [1]). From the distribu- 
tion of Pv values in Fig. [21 we find a typical uncertainty 
in the determination of Pv for a stellar population of 
a factor of ~1.4 (e.g., assuming a theoretical /3y = 1.1, 
an old stellar population may actually show a range of 
0.75 < < 1.5). This corresponds to an error in dust 
extinction, ctav-, of up to ~±0.37 mag, or magarcsec"^ 
in the case of surface brightness. If this stellar popu- 
lation is either extremely metal-poor or metal-rich, the 
uncertainty is a factor of 2, or aAv — ±0.75 mag. On 
the other hand, if an extremely young (old) stellar popu- 
lation is treated as an old (extremely young) stellar pop- 
ulation, the error in the estimated theoretical Pv value 
will become a factor of 5 or larger. This corresponds 
to a significant error in estimating the dust extinction, 
with aAv ^±l-'''5 mag. Since the extremely young and 
old stellar populations usually have distinct characteris- 
tics in the CMDs and color-color diagrams, it is highly 
unlikely that one would mistake a young stellar popula- 
tion for an old one. In §3.3 below, we will use [/-band 
observations to robustly separate pixels dominated by 
the light of younger and older stellar populations. 

3. DATA ANALYSIS 

In this section, we apply the above method to the 
observational data. Before applying it to a large sam- 
ple of galaxies, we assess its reliability in this paper for 
one galaxy: NGC 959. In future papers (Tamura et al. 
2009b, c, in preparation), we will reconstruct and analyze 
the extinction-corrected images of NGC 959 and a sample 
of 45 galaxies, spanning a wide range of elliptical-spiral 
galaxies, with GALEX FUV and NUV, optical UBVR, 
and Spitzer /IRAC images using the present method. Un- 
less indicated otherwise, the color composite images and 



the pixel coordinate maps shown in the following sections 
are oriented such that North is up and East is to the 
left. The magnitudes and colors used in th i s analysis are 
on th e AB-magnitude system (jOkd 119741 : lOke fc Giirml 
[198^ . 

3.1. Data Sets 

The galaxy we will use for our proof of concept, 
NGC 959, is a nearby late-type spiral galaxy classified 
as Sdm in the Third Reference Catalo g of Bright Galax- 
ies (RC3; Ide V aucouleurs et al. 1991"). at a distance of 
9.9 ±0.7Mpc^ ( Mould et al. 2000), with an inclination of 
50° (|Esipov et al. 1991). It has been observed from UV 
{GALEX FUV) through MIR {Spitzer /IRAC 8.0 /zm) 
wavelengths. Optical UBVR images, observed with the 
Vatican Advanced Technology Telescope (VATT), were 
taken from Taylor ct al. (2005;). The pipeline-processed 
Spitzer /IRAC images were retrieved from the Spitzer 
Archive^ via Leopard. The GALEX FUV and NUV im- 
ages were obtained from the Multi-Mission Archive at 
Space Telescope Institute^" (MAST). 

In this stud y, we us e a pi xel-based analysis, first 
introduced by iBothunI (jl986f ) in the form of pixel 
Color-Magnitude Diagrams (pCMDs). This has re- 
cently become a popular technique t o study stellar 
populations in nearby ga l axies (e.g.. [Abrah am et al.l 
[1999; Eskridg e et all 120031 : iLanvon-Foster et al. 200l 
[Welikala et al.lf2'008f) . This technique is performed in the 
same way as regular aperture photometry, but simply 
with an individual pixel as a source of flux. The advan- 
tage of using a pixel-based analysis is that this technique 
allows a two-dimensional analysis throughout a galaxy, 
or any part of a galaxy, without any object overlap or 
gaps as would be created in regular aperture photometry 
(see f igures by, e.g., IScoville et al.l 120011 : ICalzetti et all 
I2OO5II2OOI . 

To perform this pixel-based analysis on images from 
different telescopes and instruments, we first have to re- 
sample the pixel-scales and convolve all the images to 
a matching resolution. IDL^^ and IRAF^^ routines were 
used to match the pixel-scales and point spread functions 
(PSFs) of all images to the 1'.'5 pixel" ^ and ~ 5'.'3 FWHM 
of the GALEX NUV image, because these have the 
coarsest pixel-scale and PSF among the filters consid- 
ered for further panchromatic (FUV-MIR) studies (see 
Fig. [1] and Tamura et al. 2009b, c, in preparation). At the 
distance of NGC 959 (D = 9.9 ± 0.7Mpc), each 1'.'5 pixel 
spans 72 pc. As a result, the light from different stellar 
populations is mixed together within a single pixel, and 

* This value is based on its recession velocity including the 
influence of the Virgo cluster, the Great Attractor, and the Shap- 
ley supercluster, and taken from the NASA/IPAC Extragalactic 
Database (NED) as of September 24, 2008 

^ SSC: Data Archives/Analysis: 
http : //ssc . spitzer . caltech. edu/archamaly/ 

Galaxy Evolution Explorer, GR4/GR5 Data Release; 
http : //galex . stscl . edu/GR4 

IDL is distributed by ITT Visual Information Solutions (Re- 
search System Inc.), Boulder, Colorado: 
http : //rsinc . com/ldl/ 

^'^ IRAF is distributed by National Optical Astronomy Observa- 
tory (NOAO), which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation (NSF): 
http : //iraf . net 
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the observed flux ratios cannot be compared directly to 
the SSP-derived (3v values in Fig. [2l While the pixels 
with most of the light coming from older stellar popula- 
tions are not affected as much, the pixels with light from 
younger stellar populations will be affected significantly. 
Even though the light is smoothed over a much larger 
area than the area covered by a single pixel, the peak 
of the light distribution stays at the same pixel coordi- 
nate as before the smoothing. In the following, we there- 
fore statistically analyze the observed images to estimate 
the intrinsic, dust-free flux ratios for both younger and 
older stellar populations, instead of a direct comparison 
to SSP models in the previous section. We will only use 
pixels with signal-to-noise ratio of S/N > S/Nmin = 3.0 in 
all FUV-MIR filters used. 

3.2. Assumptions about the Dust Distribution 

The dust is not uniformly distributed across an entire 
galaxy. Instead, it is distributed in complicated patterns 
of wisps, lanes, and bands of thin and filamentary struc- 
tures, as well as in small clumps (e.g.. lWaller et al.lll992l : 
Deo et al. '2006|^ Using stellar radiative transfer models, 
Elmegrecn (1983) andBVitt et al. (1992) showed that de- 
pending on different spatial distributions of the dust, e.g., 
a cloud or a slab, the effect of dust extinction — includ- 
ing both absorption and scattering — varies alon g dif- 
ferent lines-of-sight. iWalterbos fc KennicuttI ()1988D esti- 
mated the variations in extinctio n values calculat e d from 
the " s ymmetry argument" (e.g., iLin dblad" 19411: lElviud 
Il956( ). ICalzetti et al.l ()1994n compared models of five dif- 
ferent dust distributions to the observational data: (1) a 
uniform dust screen; (2) a clumpy dust screen; (3) a uni- 
form scattering slab; (4) a clumpy scattering slab; and 
(5) an internal dust model (see their Fig. 8). The im- 
portance of these models is that the resulting extinction 
is different for each mo del although th e amo unt of dust 
is same. Even though iCalzetti et al.l ()1994D could not 
find satisfactory agreement between the observed data 
and these models using Large Magellanic Cloud (LMC) 
or Milky Way (MW) dust-extinction curves, models 2 
and 4 with clumpy dust distributions show in general a 
better fit to the data than the other models. We will 
assume that the dust affecting the observed light is dis- 
tributed in wisps and clumps, and is mostly in front of 
the observed stellar populations — with a much lar ger 
effect of absorption than light scattering (Bvun"1992') — 
foll owing the observ ations of, e.g.. Waller et al. (1992) 
and iDeo et all (|2006l ) . The true distribution of the dust 
would be more complex, but representable by a combi- 
nation of different types of simple geometrical distribu- 
tions, with the dust distributed not only in front, but 
also intermixed with stellar populations. Dust located 
mostly behind the stars would not be detected at the 
wavelengths considered here. Decoding the distribution 
of dust in different geometries is beyond the scope of the 
present paper, and hence deferred to a subsequent anal- 
ysis (Tamura et al. 2009c, in preparation), where we will 
have a larger sample of galaxy types and inclinations. 

Another important issue is the filling factor of dust 
within a pixel. While each pixel subtends I'.'S, the dust 
features can span either a larger or smaller area. Given a 
single average value of the dust extinction in each pixel, 
there are two extreme possibilities for the dust distri- 
bution within that pixel: (1) an extended distribution 



of a relatively thin layer of dust; and (2) small high- 
density clumps of dust covering only a fraction of the 
area in that pixel. Since the light observed in a single 
pixel is a mixture of light from different stellar popula- 
tions, the effect that we observe is a weighted average of 
the light from those stellar populations. For the former 
type of dust distribution, the dust is affecting the light 
of all stellar populations contributing to a single pixel 
equally. For the latter case, the dust is only affecting a 
small fraction of the total light in that pixel. While the 
light from stellar populations behind a dust clump is re- 
duced, the light from unextincted stellar populations is 
observed at its full strength, reducing the average effect 
of dust extinction within a pixel. If the intrinsic V-band 
flux is estimated from the observed NIR-fiux — which is 
still affected by dust extinction to some degree — this 
partial coverage of any dust extinction might cause a 
large uncertainty. Our method, on the other hand, esti- 
mates the intrinsic T^-band flux from the observed MIR 
(3.6 /im) flux, which is usually considered as extinction 
free (e.g. jFazio etal|[200l Iw'illner et al.|[200l . The ob- 
served 3.6 /im flux is therefore the total light from all 
stellar populations along the line-of-sight, ranging from 
the front to the far side of a galaxy. Since the intrin- 
sic V-band flux is estimated based on this 3.6 /zm flux, 
our method estimates the total amount of missing V- 
band flux. Hence, even though the "exact" effect of dust 
extinction depends on the true dust geometry, the to- 
tal missing y-band flux and the corresponding dust ex- 
tinction should not be affected by the dust geometry, 
unless the extinction in individual knots or filaments is 
^ 1.0 mag, thereby also affecting the observed MIR fiux, 
and preferentially affecting only young and highly con- 
centrated (on scales ^ 72 pc) stellar populations. In ret- 
rospect, this does not appear to be the case (see e.g.. 
Fig. [TT|) . although one has to be careful that this does 
not become a circular argument. 

3.3. Separating Younger and Older Stellar Populations 

Even though the effect of light-blending is significant 
in NGC 959, some pixels are still dominated by the light 
from younger stellar populations. Therefore, the first 
step in our data analysis is to separate these pixels 
from pixels whose light mostly comes from older stel- 
lar populations. Fig. [3] shows theoretical tracks of 
vs. (U-S.Gfim) using 52 Myr and 10 Gyr SSP SED 
models with a metallicity of Z = 0.008. These metallic- 
ity and ages are selected because: (1) Z = 0.008 is the 
central metallicity among five metallicities available for 
the SED library by AF03; (2) t = 10 Gyr represents an 
old stellar population; and (3) t = 52 Myr represents a 
young, but not an extremely young (t < 10 Myr) stel- 
lar population (see Fig. [2]). Since the {/-band light is 
sensitive to younger stellar populations and 3.6 /im light 
traces older stellar populations, the ({7—3. 6 /xm) color 
indicates the luminosity-weighted average age of stellar 
populations within a given pixel. The solid curve shows 
a track of the surface brightness and color as the frac- 
tion of light from the young stellar population increases. 
The fraction indicated in the figure is the mass fraction 
of stellar populations, i.e., "100%" indicates that the 
mass ratio between young and old stellar populations is 
one-to-one. A reddening vector corresponding to a vi- 
sual extinction of Ay = 1.0 magarcsec"^ is drawn in the 
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Fig. 3. — The theoretical fiy vs. ((7—3.6 fim) color-magnitude di- 
agram. The solid curve is the theoretical track produced by adding 
different fractions of a young (52 Myr) SSP SED to an old (10 Gyr) 
SSP SED: the "100%" indicates a mass-ratio between young and 
old SSPs of 1, and "200%" means the mass-ratio is two-to-one. The 
effect of an average visual extinction of 1.0 magarcsec"-^ is shown 
as a reddening vector in the lower left, and by dashed lines for se- 
lected data-points. The dotted curve traces a theoretical track for 
the color and surface brightness if the entire pixel (with a filling fac- 
tor of 100%) is affected by this reddening. The vertical dot-dashed 
line at ((7—3.6 /im) = 0.76 mag represents the color where the mass- 
weighted SEDs of young and old SSPs are contributing equally, and 
are both affected by a dust extinction of Ay = 1.0 magarcsec"^. 

lower left of Fig. [31 and is applied to selected data-points. 
Even though this is an extreme case, Fig. [3] shows that, 
once the young stellar population dominates, the color 
of mixed stellar populations tend to an asymptotic value 
of (J7— 3.6/im) ~ — 1.0 mag. In a real situation, younger 
stellar populations are known to associate with a larger 
amount of dust f e.g. Ivan HoutenI 11961: iKnapen et all 



119911: iRegan et all 120041: iBarmbv et al.M2006[ ). and the 
observed color will therefore most likely shift away from 
the no-extinction track to redder colors. 

Fig. UK shows a pCMD of /iy vs. (/7— 3.6 /im) of the 
observed images at GALEX resolution. The photomet- 
ric uncertainties plotted along the right side of the figure 
are calculated using data-points within a horizontal slice 
of fjLv ±0.1 magarcsec"^ at each point. A reddening vec- 
tor corresponding to Ay = 0.5 magarcsec"^ and assum- 
ing the LMC extinction curve is drawn at the lower left. 
The stellar populations in NGC 959 are unresolved due 
to a combination of its distance, the coarse pixel-scale, 
and a large PSF. As a result. Fig. |3K does not show 
clear separations among different stellar populations, as 
are seen for pa rtially resolved stellar populations (e.g., 
lEskridge et al.1 ,2003). Nevertheless, some branches and 
features in this pCMD are still recognizable. For sim- 
plicity, in the following sections, we will refer to pixels in 
which the flux is dominated by light from younger stellar 
populations as "younger pixels" , and pixels dominated 
by older stellar populations as "older pixels" . Visual 



inspection of Fig. 3^ (in particular, different slopes for 
groups of branching pixels on left and right sides and 
groups of pixels forming a "shoulder-like" distributions 
at — 21 magarcsec"^ in the pCMD) suggests an em- 
pirical division between the "younger" and "older" pix- 
els at (J7— 3.6/im) = 0.72 mag — close to the color indi- 
cated by a dash-dotted line at (C/— 3.6 /im) = 0.76 mag in 
Fig.[3l where young stellar populations start dominating 
the mass-fraction of a mixed stellar population within 
a resolution element and suffering a total extinction of 
Ay = 1.0 magarcsec"^. 

Before we proceed, we first perform several visual 
checks to ensure this separation of younger and older 
pixels is indeed acceptable. At colors bluer than 
(/7— 3.6/xm) = 0.72 mag, groups of pixels — or "branch- 
like features" — seem to have similar slopes as the 
reddening vector, suggesting that they represent young 
stellar populations affected by dust. This agrees with 
the notion (e.g., Ivan Houtc nl fl96lH K nap en etallllggl 
IRegan et all l2004r iBarmbv et al.l |2006D that dust is 
strongly associated with active and recent SF-regions. 
On the redder side of the pCMD, the brighter pixels 
(/iy < 21.9 magarcsec"^) form a distribution with a pos- 
itive slope (i.e., perpendicular to the reddening vector), 
indicating that a mechanism other than dust extinction 
might be affecting the fluxes in these pixels. Since the 
younger stellar populations tend to be much brighter 
than the older stellar populations (see Fig. [1]), for the 
lower surface brightness {fJ.v) pixels, the fraction of light 
from younger stellar populations as well as its corre- 
sponding dust extinction would be smaller. The theo- 
retical tracks in Fig. |2| suggest this would result in a pile- 
up of points along a fairly vertical line, e.g., the dash- 
dotted hue at (FUV — 3.6 ^m) = 0.76 mag. We therefore 
conclude that our separation of younger and older pix- 
els is reasonable, at least to first order. We will per- 
form several more checks to confirm that this separa- 
tion is appropriate (see below). At this point and at 
the present resolution, we do not see any special fea- 
tures indicating possible pixels dominated by the light 
from "extremely" young stellar populations (i.e., t<10- 
100 Myr, see Fig. [2]). Considering the effect of light- 
blending, the light from extremely young stellar popula- 
tions would most likely be diluted, and hence the separa- 
tion of pixels simply into younger and older pixels seems 
to be sufficient for this galaxy (and galaxies viewed at 
similar linear resolution). 

Fig. HId shows the spatial distribution of the younger 
and older pixels selected in Fig. [3^. The Spitzer /IRAC 
8.0 /im PAH e mission — a n ind i cator of nearby SF- 
activity (e.g., iHelou et all 12004 iCalzetti et all 120051 . 
I2007D — is overplotted in dotted contours. The distribu- 
tions of younger pixels and 8.0 /im PAH emission do not 
exactly overlap, but are shifted slightly with respect to 
one another. Since the directions of these shifts are not 
constant, they are not the result of astrometric error in 
the World Coordinate System (W CS) of the images, bu t 
are genuine features, as shown bv ICalzetti et al.r(|2005f ). 
Analysis of the PAH emission and its distribution rel- 
ative to SF-regions is beyond the scope of this paper, 
and will be deferred to a future paper. The most im- 
portant result obtained from Fig. HJa is that — while 
no coordinate information is used to select the younger 
pixels in Fig. [4^ — the selected pixels are grouped to- 
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Fig. 4. — a) A pCMD of fiy vs. (U — 3.6 fim) for observed images. The pixels dominated by the light from younger stellar populations 
("younger pixels") with ((7 — 3.6 ^m) < 0.72 mag are plotted in different colors according to the branch-like feature (see text) that they 
appear to be part of. The average uncertainty as a function of magnitude is indicated by the error bars along the right side of this panel. 
A reddening vector with Ay = 0.5 magarcsec"^ is plotted at the lower left, b) Spatial distribution of the color-coded and light gray pixels 
selected in (a). Dotted contours trace the Spitzer 8.0 fim PAH emission. The younger pixels form contiguous regions; they are not randomly 
distributed. Different branch-like features in panel (a) correspond to regions at both systematically different galactocentric radius and with 
systematically different PAH emission. The younger pixels on the red-coded branch correspond to regions with strong PAH emission, and 
likely significant Ay, that are also relatively close-in; the purple pixels are faint at 8.0 fim, suffer little extinction, and are found at larger 
distances from the center. 



gether into cohesive regions that coincide, or are close to 
peaks in the IRAC 8.0 /im emission. Moreover, different 
branch-like features in Fig. 3^ turn out to correspond to 
regions at both systematically different distances from 
the galaxy center and systematically different PAH sur- 
face brightness. The younger pixels that form the high- 
surface brightness (fiv) feature that is color-coded red 
in Fig. [4^ originate mostly from two regions that are lo- 
cated near one of the strongest peaks in the 8.0 fim emis- 
sion and likely suffer significant extinction. Their high 
surface brightness appears mostly due to the relatively 
small distance from the galaxy center and exponential de- 
cline in surface brightness of the disk of NGC 959 (e.g., 
iHeraudeau fc Simienlll996l : iTavlor et al.i r2005). Younger 
pixels on progressively lower surface brightness features 
(color-coded in order: orange, gold, green, blue and pur- 
ple) correspond to regions with progressively larger dis- 
tance from the center and fainter PAH emission (smaller 
Av). 

To further confirm that our selection of the younger 
pixels is not a random result, we also visually compare 
Fig. [4|3 to color composite images of NGC 959. Fig. [5] 
shows UVR color composites at two different spatial res- 
olutions. Both are composed from the same U -, V-, 
and fl -band images observed at the VATT (Tavlor et al.l 
I2005D . Fig. is presented at the original pixel-scale 
of 0'.'37 pixel^^, with the PSFs in all images matched 
to - 1'.'3 FWHM. In Fig. [Hd, the pixel-scale and resolu- 
tion were matched to that of the GALEX NUV image, 
i.e., 1'.'5 pixel-i and 5'.'3 FWHM. Spitzer /l^kC 8.0 /im 
contours (green) are overlaid in the latter image for a 



comparison to the pixel-map (Fig. |4|d). Wc find that the 
spatial distribution of the selected younger pixels clearly 
follows that of bluer regions in NGC 959. 

3.4. Selecting the Theoretical V -to-3.6 fim Flux Ratio 

Having separated the "younger" and "older" pixels, 
we now want to estimate the theoretical extinction-free 
flux ratio for each pixel. As mentioned above, however, 
the theoretical (3v value cannot be simply selected from 
Fig. [21 due to the smoothing over stellar populations and 
subsequent blending of light. For pixels dominated by 
the light from older stellar populations, this should not 
be a major problem, since their theoretical flux ratios 
do not change much with age (see Fig. [2]). The prob- 
lem is the effect on pixels dominated by younger stellar 
populations. While a young SSP has a theoretical flux 
ratio of ^ 5 (or up to ~ 11 for an extremely young stellar 
population), the extinction-free (Sy value for mixed stel- 
lar populations depends strongly on how large a fraction 
of the light originates from underlying and neighboring 
older stellar populations. Since we cannot measure the 
exact fractions of light from younger and older stellar 
populations, we have to find another way to estimate 
theoretical (3v values. 

We first calculate the flux ratio of the observed V 
and 3.6 /im images for each pixel to characterize the ef- 
fects of hght-blending. The result is plotted in Fig. [6l 
where the white-colored histogram represents the distri- 
bution of f3v values for the pixels dominated by the light 
from older stellar populations ("older" pixels), and the 
gray-colored histogram for those dominated by the light 
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Fig. 5. — UVR color composite images of NGC 959 at two different spatial resolutions, a) A color composite image constructed at the 
original pixel-scale of the VATT CCD, 0^37 pixel" ^ with a seeing of I'.'S FWHM. Dust lanes and stellar populations with different color 
are readily discernable in this image. 6) A color composite image constructed from images that are resampled to the GALEX pixel-scale 
of 1'.'5 pixel"^, and convolved to the GALEX PSF of ~5''3 FWHM. IRAC 8.0/im emission is indicated by green contours. Compared to 
the higher resolution image, dust lanes and stellar populations in this lower resolution image are not as evident. However, the blue young 
stellar populations remain clearly visible in the lower resolution image. 



from younger stellar populations ("younger" pixels). The 
older pixels have a peak around /3v' — 1-0 with a sec- 
ondary peak, or a shoulder, at j3v — 1-25 . This confirms 
our assumption that light-blending has only a minor ef- 
fect on older pixels. The younger pixels display a peak at 
(3y ~ 1.32, which is much lower than the theoretical ratio 
(4 < /3y < 7, see Fig. [2|) . This indicates that younger pix- 
els are significantly affected by the light from underlying 
and neighboring older stellar populations. An important 
feature of the distribution for the younger pixels is that 
the observed range of j3v values is narrow and concen- 
trated within 0.9 < /Jy < 1.6 . This implies the effect of 
light-blending is rather uniform and consistent for these 
pixels. Stated differently, the contaminating older stel- 
lar population is distributed much more smoothly than 
the younger stellar population (see Fig. [51 Fig. [T3] and 
Figs.IUHIZl). 

Another important feature of Fig. [H] is that the tail 
of the distribution toward lower /3y values is larger than 
the tail toward higher values. Since dust extinction is the 
primary cause for the reduction of the l/-band flux, the 
pixels with anomalously low /3y values are most likely 
located in the high-extinction regions within the galaxy. 
While both groups of pixels have low-/3y tails, the rel- 
ative size of the tail, compared to the size of the main 
distribution, is much larger for the younger pixels. This 
indicates that dust extinction is more significant for the 
younger pixels. 

Because interstellar dus t is usually conc entrated in rel- 
atively small regions (e.g.. lDeo et"aLll2006f) — while most 
other regions suffer minimal dust extinction — we should 
be able to estimate the intrinsic dust-free V^-to-3.6/im 
flux ratio, /3v,0i for both stellar populations from Fig. [51 
While the peak of the distribution for the older pixels 
is at (iv = 1.04, this value occurs toward the lower end 
of the distribution. A statistical analysis with 2a clip- 



ping shows that the mean is at (3vfl = 1.10 with a stan- 
dard deviation of 0.14. Even though the exact value of 
the intrinsic flux ratio varies from pixel-to-pixel due to 
the different stellar population components in each pixel, 
we adopt this value of fiv^ = 1.10 ± 0.14 as the dust-free 
flux ratio for the older pixels. Unlike the older pixels, 
the pv distribution of the younger pixels is more asym- 
metric, with a much larger tail toward lower (3v val- 
ues. After 1(7 clipping, the statistical mean is pv = 1.29, 
which is smaller than the peak value of /3y = 1.32 . Be- 
cause younger stellar populations are more affected by 
dust than older stellar populations, which is also appar- 
ent from the much larger and wider lower-/3y tail, we 
elected to use /3v,o_ob = 1-32 t^^s (the peak /?y value) 
as the extinction- free flux ratio for the younger pixels. 
The added subscript, OB, indicates that the value is for 
the pixels dominated by the light from younger stellar 
populations, but does not imply that such pixels have no 
light contributed by older stellar populations. 

3.5. Calculating the Flux Difference 

Having determined appropriate extinction-free flux ra- 
tios, the next step is to calculate the y-band flux dif- 
ference, A/y, for each pixel. Using the theoretical flux 
defined in equation ([2]), the flux difference is calculated 
as: 

A/y-[/y-(/3yoX/3.6^m)], (4) 

where, fv and /3.6/jm are the observed pixel-fiuxes in the 
optical V and IRAC 3.6 ^m bands, and /3y,o is the esti- 
mated extinction-free fiux ratio. The distribution of A/y 
values is shown in Fig. [3 for both younger and older pix- 
els. The solid histograms are for the A/y values calcu- 
lated from our adopted values for /3yo, while the dashed 
(dotted) histograms are computed for /3yo values at the 
upper (lower) bound of the quoted la uncertainty range 
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Fig. 6. — Histograms of observed V-to-3.6 fim flux ratio, /By, for 
individual pixels dominated by the light from older stellar popu- 
lations ("older pixels") and younger stellar populations ("younger 
pixels"), as selected in Fig.|3] For both populations, the low-ratio 
tail of the distribution — presumably due to dust extinction — is 
more extended than the tail on the high-ratio side. We adopt an 
extinction- free flux ratio, PvO) of 1.10 ±0.14 for the older pixels 
and of 1.32 ^ q ?5 for the younger pixels. Once 2a clipped mean and 
peak values are selected for older and younger pixels, respectively, 
corresponding uncertainties are determined from the FWHM of 
iSy distributions. Because of light-blending from different popula- 
tions contributing to the flux in a given pixel, the /By values for 
the young stellar population are significantly reduced compared to 
those calculated for SSP models in Fig. [2] 

in Fig. [HI While Fig. [H] shows simply the level of the 
l/-band flux relative to the flux in 3.6 /im fllter, Fig. [7] 
shows the absolute difference between the observed and 
estimated theoretical V^-band pixel- fluxes. Regardless of 
the selection of Pv,o, both histograms — for both younger 
and older pixels — show that the number of pixels with 
positive Afv goes to zero quickly, and that there are 
distinct tails extending to large negative values of A/y. 
For the younger pixels, the effect of the uncertainty in 
the selection of Pv.o is small (i.e., the three histograms 
show similar distribution in Fig. Uh)- The older pixels, 
however, are affected more (Fig. [7^). Since older stel- 
lar populations tend to have less dust intermixed, and 
hence suffer less extinction, the dashed histogram for the 
larger value of Pv^ produce an excess of pixels with large 
values of Ay- On the other hand, the dotted histogram 
for the smaller value of Py^ will produce many pixels 
with an unphysical excess of visual flux, indicating that 
this value must be a robust lower bound on /3v,o- This 
confirms that /3y.o = 1.10 is a reasonable and appropriate 
value for the older pixels. 

Next, we check the spatial distribution of the calcu- 
lated A/y values. If the distribution of negative Afy 
follows genuine galactic features, such as SF-regions, spi- 
ral arms, and PAH emission, then this strengthens our 
argument that our method largely traces the dust ex- 
tinction. The spatial distribution of Afy is shown as a 



pixel-coordinate map in Fig. [SI with the IRAC 8.0 /Ltm 
emission overlaid as dotted contours. This map demon- 
strates that the distribution of Afy values is not ran- 
dom at all, but closely associated with genuine galactic 
structures. Visual comparison of this map to the color 
composite image in Fig. [S)d confirms that pixels with a 
large deficiency of observed F-band flux are not simply 
corresponding to pixels appearing darker (lower surface 
brightness). Instead, pixels with Afy < are distributed 
around the regions that appear bluer, as well as near 
the center of the galaxy. The higher resolution image 
of Fig. [S^ shows that some pixels with larger negative 
values of Afy (darker gray in Fig. [5]) trace the visible 
dust lanes, seen in silhouette, and the bluer SF-regions. 
This demonstrates that even though some of the dust is 
not visually conspicuous (as in the lower resolution im- 
age in Fig. [5)3), our method is capable of estimating dust 
extinction and its spatial distribution from observations 
in only two broadband filters — with a third filter ([/- 
band) serving only to robustly separate younger pixels 
from older pixels (Figs. [3] and [4]). Since regions with vis- 
ible or plausible dust content are recovered well by our 
method, this gives us confidence that regions for which 
this method indicates a low dust content could also be 
real. 

Before we proceed to estimate the amount of dust 
extinction in each pixel, there is one more check that 
we must perform to support our separate treatment of 
the younger pixels. The main panel of Fig. [9] shows 
a map of Afy for the younger pixels for our adopted 
value of /3y,o_OB = 1-32, while the inset shows the re- 
sult when these pixels are treated as older pixels with 
Pv,0-OB — Pv,o = 1.10 . The lighter shade of gray of the 
pixels in the inset panel indicates that most of them 
have positive Afy, suggesting that the observed F-band 
flux is equal to or larger than expected from the 3.6 fim 
data, i.e., that the younger pixels miss either no V- 
band flux or show an unphysical excess of observed V- 
band flux. Since young stellar populations are usually 
associated with relatively large amounts of dust (e.g. , 
van Ho utcn 19'61: Knapen et al.|[T99]] : IRegan et al.ll2004l : 
Barmbv et al. 2006 ). a value of Py_o = 1-10 is clearly too 
small for the younger pixels. Our separate treatment 
of the younger pixels, adopting /3y,o_OB = 1-32 to'is as 
outlined above, is therefore appropriate. For that larger 
value of /3v,o, the younger pixels are found to suffer small- 
to-moderate amounts of extinction by dust, consistent 
with the known presence of dust at or near these regions 
of recent star formation and producing consistent results 
with previous studies. 

3.6. Measuring the Dust Extinction 

We can now estimate the most likely amount of 
dust extinction in each pixel. With /3y^o — 1.10 
(and /3v.o_OB = 1-32) and AB-magnitude zeropoint, 
^P,AB = 48.59 mag, equation ^ becomes: 

Ay = my - [-2.51og(/3v,0 X /s.GMm) - V^p^As] ■ (5) 

We cannot simply apply equation ^ to all pixels, how- 
ever, since some pixels have (3y > Py^ (see Fig. [5]) and 
Afy > (see Fig. [7]), i.e., seemingly implying an unphys- 
ical negative dust extinction. Ay . Because the f3y^o value 
is the estimated extinction-free flux ratio, rather than the 
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Fig. 7. — Histograms of Afv, the difference between the observed and theoretical V-band flux for selected pixels with Ught dominated 
by (a) older and (6) younger stellar populations. In each panel, the solid histogram is for the adopted value of /3v,0) while the dashed 
(dotted) histogram assumes a value for /3v,o that corresponds to the upper (lower) bound of the uncertainty range determined in Fig. |6] 
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Fig. 8. — Pixel map of the flux difference, A/\a, rendered at 
the resampled pixel-scale of I'.'S pixel" ^. Darker gray-scale levels 
indicate a larger deficiency of observed V-band flux compared to 
the modeled intrinsic flux. Black dotted contours trace the IRAC 
8.0 fim PAH emission. While it is hard to visually trace the distri- 
bution of dust in Fig. the locations of the dust lanes in Fig. [SJi 
generally agree with the distribution of the flux-deficient pixels in 
the present figure. 

true dust-free ratio, it is possible that some pixels have 
an observed /3y value that is larger than Pv,o- We there- 
fore have to be careful how we treat these pixels in our 
analysis. 
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Fig. 9. — Comparison of the flux difference, Afv, inferred for 
pixels dominated by younger stellar populations. The main panel 
shows the result when the younger pixels are treated separately, 
adopting /3v,o_OB = 1-32 . Here, most pixels show that they are 
missing V-band flux, as expected. The inset panel shows the result 
for no such special treatment (i.e., assuming f3vo_OB = f^vo = 1.10, 
as for pixels dominated by old stellar populations). Here, most of 
the younger pixels show an unphysical excess of observed V-band 
flux. Black dotted contours trace 8.0 fim PAH emission. 

As mentioned above, the pixels with positive A/y 
(Fig. [8]) appear darker than neighboring regions in the 
color composites of Fig. [5] Considering also that the 
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8.0 fim PAH emission is weak in these regions (i.e., these 
pixels are surrounded by outer 8.0 /xm contours with 
large spacing between neighboring contours), these pos- 
itive Afv pixels are not caused by SF-activity, but 
are most likely the result of underestimating the in- 
trinsic f3v,o value. The actual amount of dust extinc- 
tion might also be minimal or zero {Ay — 0) for these 
pixels. We adopted a single Pv,o for younger and 
for older pixels, yet the intrinsic flux ratio for each 
pixel will vary slightly based on different factors. As 
a result, some pixels with positive A/y are expected, 
even in the absence of noise. For our adopted Pv,o of 
1.10 (1.32) for older (younger) pixels, most pixels with 
excess flux have A/v < 1.0 x 10^^^ crgscm^^s^^, and 
only a very small fraction has an excess as large as 
~2.0x 10~^^ ergscm~^s~^ (see Fig. [7]). In the follow- 
ing, we will therefore assume that these pixels suffer no 
extinction, i.e., that Ay =0 magarcsec"^. 

3.6.1. Impact of the Uncertainty in l3v,o 

Since neither the relationship between Ay and Afv 
nor the relationship between Afy and py is linear, the 
la errors in the estimated extinction-free flux ratios 
i'^l^v.o ^ ±0.14 and (Ti3v,o-ob = +0.09 & —0.15) cannot be 
simply converted to corresponding values. Instead, 
we will assess how the distribution and values of Ay 
change, as we vary the estimated dust-free flux ratio from 
the adopted value of Py.o to Py.o ± ct/3^„. 

Fig. [TO] shows the spatial distribution of the esti- 
mated dust extinction. Ay, for each pixel in NGC959, 
for two sets of theoretical extinction-free flux ratios. 
Fig. [TOk shows the distribution of Ay inferred for our 
adopted Py^o values of 1.10 and 1.32 for older and 
younger pixels, while Fig. [TOb shows the result for Py^ 
values set at the higher end of the uncertainty range 
(i.e., /3v_o = l-24 and 1.41, respectively). The mean 
and maximum dust extinction for all pixels in Fig. llOb 
(Fig.[lOlD) are Ay,„,ean = 0.064 (0.15) and Ay,„,ax = 0.80 
(0.93) magarcsec"^. As expected, the dust-extinction 
map (not shown) for the lower limit, Py.o—f^Pv oy is cov- 
ered mostly by Ay^O pixels. These results are sum- 
marized in Table 1, where the last column indicates the 
fraction of pixels with Ay = mag arcsec"^ in the image. 

Even though the estimated Py^ value changes, Fig. [TO] 
shows that the distribution of dust extinction follows 
the structures of the galaxy and the 8.0 /im contours. 
Fig. [TO] and Table 1 also show that the effect of vary- 
ing the theoretical Py^o value is not equal to simply 
adding or subtracting a constant value AAy to the dust- 
extinction values calculated for our adopted value of 
Pyo- As Pyo change from 1.10 and 1.32 (Fig. [TOJi) 
to 1.24 and 1.41 (Fig. [TOb). some pixels near pixels 
with Ay > magarcsec"^, which originally were deemed 
extinction-free, now suffer a slight amount of dust extinc- 
tion. Other pixels that are further away, e.g., pixels in ar- 
eas that are faint at 8.0 /zm, stay at Ay = mag arcsec"^. 
This confirms our assumption that these pixels have min- 
imal or no dust extinction. 

3.6.2. Analysis at Higher Spatial-Resolution as 
Confirmation 

While individual dust features usually are small-scale 
structures, our initial analysis was performed at the low 



spatial resolution of the GALEX NUV image. To deter- 
mine if our result is a true measure of dust extinction 
or not, we repeat the same analysis with higher spatial 
resolution images, using only U , V , and IRAC 3.6 /im 
images. Since the 3.6 /im image has the coarsest pixel 
scale of T.'2 pixel"! ^ pgp ^j^j^ _2'.'2 FWHM, the 
ground-based U- and F-band images are registered, re- 
sampled, and convolved to match the orientation, pixel 
scale, and resolution of the 3.6 /im image. Fig. [TT] shows 
the distribution of estimated Ay values at this higher 
spatial resolution. The 8.0 /im contours are overplotted 
at the native IRAC resolution of - 2'.'3 FWHM. 

An important difference between the two spatial res- 
olutions is that, while the maximum dust extinction is 
^y,max=i0.8 magarcsec"^ in Fig.[TO]i, Ay^^ax in Fig.[TT] 
reaches ~2.3 magarcsec"^. This jump in Ay value 
is expected, since the coarser GALEX PSF smoothes 
out the effect of dust extinction and reduces the av- 
eraged Ay for each pixel. Since the ratio of the ef- 
fective areas of the GALEX NUV and Spitzer 3.6 /im 
PSFs is (5'.'3/2'.'2)2c:i6, the typical extinction per pixel 
should change, to first order, by ^^1.9 magarcsec"^. 
Other factors, such as uncertainties in estimating Py^Q 
and blending of light from structures that are unre- 
solved with the coarser while resolved with the smaller 
PSF, also play a role. The observed shift in Ay of 
AAy ~ 1.5 magarcsec"^ is therefore broadly consistent, 
while the 0.4 magarcsec"^ difference illustrates the ef- 
fects of the dumpiness of stars and dust on scales of 
2'.'2-5'.'3 (~ 110-250 pc at the distance of NGC 959). The 
higher spatial resolution images are better at tracing 
smaller dust features and their higher extinction values 
within a single pixel. Nonetheless, the overall distribu- 
tion of dust extinction is similar in Figs. [TO] and [TT] While 
the fine details of the measurable dust extinction and its 
spatial distribution depend on the resolution, both ex- 
tinction maps trace the 8.0 /im PAH emission and trace 
genuine galactic structures. 

In conclusion, from the series of tests described above, 
we find that the two-dimensional distribution of dust ex- 
tinction can, indeed, be reliably estimated by our pixel- 
based method and the observed V-to-3.6 fim flux ratio. 

3.6.3. Interesting Regions 

Having produced maps of the spatial distribution of 
dust extinction in NGC 959, we now discuss some of the 
most interesting dust features in Fig. [TO] (and Fig. [TTjl . 
While the distribution of regions with larger Ay closely 
traces the SF-regions, a not previously identified bar, and 
arm-like structures, there are several regions that draw 
our attention. These are: (1) a compact region at the 
Northeast edge (NE; upper-left) of NGC 959 that appears 
to suffer very high extinction; (2) another such high- 
extinction region at the Northwest edge (NW; upper- 
right); and (3) an extended area of moderate dust extinc- 
tion along the Southeast (SE; bottom), rim of the galaxy. 
Since all of these features are also visible in Fig.[TTl these 
must be real features. To visually confirm small-scale 
dust features in NGC 959, and check their interpretation 
as genuine galactic features as opposed to chance super- 
positions of unrelated objects, we created another color 
composite from Archival HST /WFFG2 F450W, F606W, 
and F814W images, shown in Fig.[l2l The two circles in 
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Fig. 10. — Spatial distribution of Ay, inferred for Pv.o values of (a) 1.10 and 1.32 (mean), and (b) 1.24 and 1.41 (upper limit), for older and 
younger pixels, respectively. The same gray-scale is used for both (a) and (6). The pixels with the lightest shade of gray have no detectable 
extinction (Ay = mag arcsec"^). The darkest pixels have a visual extinction Ay ~0.8 mag arcsec"^ in (a) and ~0.93 mag arcsec"^ in 
(6). Dotted contours trace the 8.0 /im PAH emission. These extinction maps indicate that the higher extinction coincides with SF-regions 
that appear blue in the color composite images (see Fig. [5}, as well as with several dust features visible in Fig.[SJi. 



TABLE 1 

Dust Extinction Based on Different f3v,o Values 



Selected ^v^^^^^ ^v^^x % of pixels 

/3v,o Level (^v,o /3v,o_OB (mag arcsec"^) (mag arcsec"^) with Ay =0 

Adopted /3vo 1-10 1.32 0.064 0.80 45.0% 

Upper (+ o-)' 1.24 1.41 0.15 0.93 16.4% 

Lower (-cr) 0.96 1.27 0.015 0.63 83.1% 

Mean Ay value of all pixels analyzed for the galaxy. 



this image mark regions (1) and (2). 

The NE high- Ay region, region (1), is centered 
around pixel coordinates [x, y] ~ [128, 169] in Fig. [TOl 
Even though nothing conspicuous is visible in the low- 
resolution color composite (Fig. [5^), the HST image 
(Fig. [121) reveals a compact, bright red source at the cen- 
ter of the marked region. Without morphological indi- 
cators, multi-filter photometry, or spectroscopic informa- 
tion for this particular object, it is hard to decide whether 
this is a reddened stellar population within NGC959's 
disk, or a background (foreground) object that is visible 
through (against) the disk. We do note, however, that 
the color of this object is very similar to that of the edge- 
on background galaxy that is visible at the bottom left 
of Fig. [T2l Region (2) is located in the NW corner of 
the galaxy around [a;, y] ~ [178, 155]. There is no F606W 
coverage for this region in Fig. [121 but there is no object 
discernable in the other two HST filters. At present, we 
lack sufficient information to establish if these regions 
are truly associated with NGC 959. Until further evi- 
dence is obtained, we will treat these regions as a part of 
the galaxy. 

The pixels comprising feature (3), the moderate- 



extinction region running along the SE rim of the galaxy 
disk, initially did not show up as having particularly 
large-negative Afy values (see Fig. [51). Once the dust 
extinction is calculated, these pixels do stand out with 
low-moderate Ay values, indicating the presence of an 
extended dust structure. The higher-resolution dust ex- 
tinction map (Fig. fTI]) also shows that the pixels in this 
region suffer higher dust extinction than inferred for the 
opposite (NW) rim of the galaxy. Visual comparison 
to Fig. [121 confirms that a faint dust lane runs along 
the SE rim, apparently tracing an outer spiral arm or 
armlet. This shows that our dust extinction measure- 
ment using the y-to-3.6 /im flux ratio is sensitive to 
even low amounts of dust extinction per resolution el- 
ement, whether inherently smoothly distributed or re- 
sulting from a small filling factor. 

3.6.4. Correcting for Dust Extinction 

Finally, using the calculated Ay values for each pixel, 
we correct the observed l^-band image of NGC 959 to 
reveal the true underlying stellar populations. The 
observed (uncorrected) and extinction-corrected images 
are shown at GALEX resolution in the top panels of 
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Fig. 11. — Distribution of dust extinction at the (higher) spatial 
resolution of the Spitzer/IUAC 3.6 image {1'.'2 pixel-l and 
~2'.'2 FWHM). The dotted contours trace 8.0 /im PAH emission at 
a native resolution of '~2"3 FWHM. This map is presented in the 
instrument coordinate system, oriented as indicated at the upper 
left. 
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nent, while other galactic structures (particularly the 
central bulge) also become better defined. Especially 
for areas SE (centered on [x, y] ~ [142, 148]) and NW 
{[x, y] ~ [159, 155]) of the galaxy center, the F-band sur- 
face brightness becomes much brighter in areas where 
dust lanes are evident in the color composites (Fig. [S^ 
and Fig. [T^ . The bottom two panels show the observed 
Spitzer /IRAC 3.6 ^m and 4.5 /um images (Figs. 13c and 
13d), which trace the di stribution of the underlying older 
stellar populations (e.g.. lRegan et al1l2004l : rWillner et al.l 
|2004). These MIR images and the extinction- corrected 
y-band surface brightness distributi on show exc ellent 
qualitative agreement. Therefore, as iReganI ()200Q) did 
with optical-NIR images and radiative transfer model- 
ing, we successfully corrected the dust extinction with 
images in only two filters {V and 3.6 /.tm) — with a U- 
band image serving only to robustly separate pixels dom- 
inated by the flux from younger stellar populations from 
those dominated by older stellar populations. 

4. APPLICATION TO OTHER FILTERS 

Given the amount of visual dust extinction. Ay, in 
each pixel, we can calculate the extinction in any other 
filter. The extinction in a given filter depends on its 
throughput as a function of wavelength, as well as on 
the metallicity of the stellar populations of the galaxy 
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In a forthcoming paper (Tamura et al. 2009b, 
in preparation), we will use UV-MIR multi-filter dust- 
corrected surface photometry for a detailed analysis of 
the stellar populations within NGC959. 




Fig. 12.— Color composite of Archival ff5'r/WFPC2 F450W, 
F606W, and F814W images. The green cir cles indicate the two 
regions with high apparent Ay in Fig. 1101 While there is no 
object evident within the right green circle, the circle in the up- 
per left corner contains a red object. Along the southern rim of 
NGC 959, a faint dust lane is discernable that appears to trace 
an outer arm. The typical extin ction along this dust feature is 
Av~0.3-0.4 magarcsec"^ in Fig. 1101 

Fig. [131 using the same gray-scale lookup table. The 
surface brightness distribution before extinction correc- 
tion (Fig. [T5b ) has lower contrast overall and fewer high- 
contrast features. Most conspicuous after applying our 
dust correction (Fig. [13b ) is that bluer regions in the 
color composite of Fig. [Sja become much more promi- 



4.1. Extinction Curves 

In general, dust extinction increases toward shorter 
wavelengths. For metal-poor stellar populations — such 
as those in the Small Magellanic Cloud (SMC) — the 
dust-extinction curve is largely monotonic as a func- 
tion of wavelength. The shape of the extinction curve 
becomes more complicated for metal-rich stellar pop- 
ulations. Especially the extinction at shorter wave- 
lengths (A < 2200 A) may be significantly affected by 
the 2175 A interstellar dust feature and the increased 
steepness of the UV extinction (e.g.. I Calzetti et 311119941 : 
Gordon et al. 2003, see also Fig. [Til here). 

[Gordon et al.l (|2003l) compared the observations of 
many stars to calibrate the extinction curves — the wave- 
length dependent extinction relative to that at 0.55 /xm 
— for SMC, L MC, and MW type du st. Based on Table 4 
(and Fig. 10) of lGordon et al.l (|2003[ ). we interpolated the 
published extincti on curves on a finer wavelength grid. 
In the NIR, G ordon et al.l (j2003) have only one data- 
point for each of the 2MASS JHKs filters, and the SMC 
extinction curve in particular appears quite uncertain. 
To cover longer wavelengths up to the Spitzer /IRAC 
filters, in the NIR an d MIR, we adopte d the Galac- 
tic extinction curve of iFitzpatrickl (|1999[ ). rescaled for 
R = Ay /E{B — V) = 3.1 . The adopted extinction curves 
for SMC, LMC, and MW metalhcities are plotted in 
Fig. [m For wavelengths shorter than the 2175 A in- 
terstellar dust feature, the SMC and MW type extinc- 
tions differ by ~ 2 mag. The bottom panel of Fig. [TH 
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Fig. 13. — Comparison of observed (top left) and extinction-corrected (top right) V-band images and the observed Spitzer /IRAC 3.6 
and 4.5 /^m images (bottom) of NGC 959 at GALEX resolution. The MIR images are good tracers of the underlying stellar distribution. 
The extinction-corrected image in V displays more brighter (darker shade) pixels around blue SF-regions (see Fig. [SJ and better resembles 
the MIR images than does the observed V-band image. 



shows the relevant filter throughput curves for compari- 
son. The amount of dust extinction in each filter, Afiitcr, 
is calculated as a ratio to the extinction in V: 



^filter My 



ATfiite..(A) [A{\)/Av]d\ 
Xx7fiitcr(A)dA 



(6) 



where, Tfiiter(A) is the throughput curve for each filter, 
and [y4(A)/^y] denotes an extinction curve (Fig. [H]) . 
Table 2 summarizes the computed dust extinction for 
the different extinction curves. For the 2MASS JHKg 



filters, the extinction calculat ed from the SMC, L MC, 
and MW extinction curves of lGordon etld] (l200l dif- 
fers significantly (indicated with asterisks in Table 2). 
Since the extinction at longer wavelengths is progres- 
sively insensitive to metallicity, we conclude that these 
values must be hi ghly uncerta i n. In stead, we use the 
extinction curve of iFitzpatrickl ()1999D to provide an up- 
per limit to the dust extinction in the filters longward of 
1 /im (parenthesized values in Table 2). Since the extinc- 
tion is small (Afntcr/^y <0.1), following previous stud- 
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Fig. 14. — Top panel: Extinction curves for SMC (pluses), LMC 
(open diamonds), and MW (crosses) metallicity from FUV to MIR 
wavelengths. The smooth extinction c urves were interpolated from 
the data points ofl Gord on et al.l 02OO3). Also shown is the Galactic 
extinction curve from Fitzpatrick (1999), which we adopt for the 
NIR and MIR filters. Bottom panel: Total throughput curves, 
T(A), for different telescope-filter combinations (sec also Fig. [TJ. 
The GALEX filter curves were scaled up by a factor of 10 for 
better visibility. 

ies fe.g., lFazio et al.|[200^ iRegan et al.ll2004D . we assume 
that there is no measurable dust extinction in the IRAC 
filters. 

An important assumption for the extinction curves de- 
scribed above is that the dust is distributed in the form of 
a "diffuse-ISM" or in a diffuse-screen geometry, which is 
applicable only for nearby stars and star clusters within 
our own Galaxy. In extragalactic objects, the dust ap- 
pears to be d istributed in smaller clumps of much higher 
density (e.g.. lDeo et al.li2006.) intermixed with the stars. 
Even though the effect of dust extinction is averaged (or 
smoothed out) within an aperture or a resolution ele- 
ment — a single pixel in our pixel-based analysis — the 
properties of dust extinction (e.g., the 2175 A feature and 
the steepness of the UV-extinction curve) are governed 
by the actual geometry of the dust distribution within 
a galaxy. This means that two regions with the same 
average visual extinction Ay can have different amounts 
of extinction A\ at another wavelength, either because 
of differences in metallicity or because of different dust 
geometries. 

Two extreme cases of dust geometry are: (1) a uni- 
form thin slab (i.e., commonly referred to as a "diffuse 
ISM"); and (2) dense clumps covering a small fraction 
of a resolution element. While we would like to perform 
a detailed study of the dust geometry and its effect on 
the extinction curve for extragalactic objects, this is be- 
yond the scope of the present paper. Here, we briefly 
discus s the effects of the different dust geometrie s stud - 
ied bv lWhittet eFall (l200ll . [200l . IWhittet et al] (l200lh 
studied a total of 27 sight-lines (stars) toward the dark 
clouds in the Taurus region. One of the results from their 
study is that Ry = Ay / E{B — V) changes from a "nor- 
mal" to a "dense cloud" regime once the extinction ex- 
ceeds a threshold value of Ah — 3.2 mag. IWhittet et al.l 
(|2004f) subsequently studied the effect of dust geometries 
in detail, which included a thin "diffuse screen" and a dif- 



TABLE 2 
Dust Extinction (Aflnor/Av) 



Instrument 


•Reenter 


X 








& Fihcr 


(Atm) 


(Atm~i) 


SMC 


LMC 


MW 


GALEX FUV 


0.153 


6.536 


4.56 


3.16 


2.70 


GALEX NUV 


0.227 


4.405 


2.91 


2.59 


2.53 


VATT U 


0.360 


2.778 


1.65 


1.64 


1.50 


VATT B 


0.437 


2.288 


1.37 


1.35 


1.29 


VATT V 


0.542 


1.845 


1.00 


1.00 


1.00 


VATT R 


0.642 


1.558 


0.81 


0.81 


0.81 


2MASS J 


1.235 


0.810 


0.14* 


0.29* 


0.36* 












(<0.26) 


2MASS H 


1.662 


0.602 


0.15* 


0.10* 


0.18* 












(<0.17) 


2MASS 


2.159 


0.463 


0.03* 


0.09* 


0.04* 












(<o.ii) 


IRAC 3.6 Atm 


3.550 


0.282 






(<0.06) 


IRAC 4.5 Atm 


4.493 


0.223 






(<0.04) 


IRAC 5.8 Aim 


5.791 


0.173 






(<0.03) 


IRAC 8.0 Aim 


7.872 


0.127 






(<0.02) 



Note. — An asterisk indicates that the uncertainty is sig- 
nificant compared to the actual extinction values. The values 
in parenth esis are calculated using the Galactric extinction 
curve from lFitzpatrickl l ll999l) . 



fuse screen with an embedded "dense cloud" (see their 
Fig. 1). For a detailed analysis and discussion, we re- 
fer the reader to their papers. The main effect of the 
"dense cloud" geometry on the extinction curve is to 
w eaken or remove the 2 175 A extinction bump (see Fig. 5 
of IWhittet et al.ll2004f) while having little effect on the 
extinction curve at other wavelengths, which remains at 
th e same level as for the mean "diffuse ISM" (Fig. 2 
of'Whittet et al."2004l. The ex tinction curves recreated 
from I Gordon et al. (2003) and iFitzpatTick (1999), and 
the calculated extinction values Afntei/Ay (Fig. [T4l and 
Table 2) are therefore treated as the upper limit to the 
dust extinction from different dust geometries. Since the 
2175 A feature is covered only by the GALEX NUV filter, 
the uncertainty associated with different dust geometries 
is assumed to be minimal in all filters except the GALEX 
NUV. Since the GALEX NUV filter is not used in our 
method, we will defer further analysis of the NUV filter 
to subsequent papers (Tamura et al. 2009b, c, in prepa- 
ration). 

4.2. A Dust-Free View of NGC959 

Since NGC 959 is classified as an Sdm galaxy in the 
RC3, we assume that it has an average metallicity be- 
tween the SMC and the MW value. We therefore adopt 
the LMC extinction curve to estimate the extinction in 
filters other than V . Using the LMC extinction ratios 
Af^\tev/Ay from Table 2, we scale the F-band extinction 
for each pixel as: 

Alter = Ay X (AfiitcrMl^) , (7) 

where. Ay was computed using equation ([5]). Figs. lTSl lTfl 
show the uncorrected (left panels) and the extinction- 
corrected images (right panels) for NGC 959 from the 
FUV to R at GALEX NUV resolution. In all fihers, the 
SF-regions become more clearly visible in the extinction 
corrected images. Compared to the optical l^-band and 
the MIR 3.6 and 4.5 /im images (Fig. IT5| . in which the 
highest surface brightness is reached in the galaxy center, 
the brightest regions in the extinction-corrected GALEX 
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Fig. 15. 



Comparison of observed (left) and extinction-corrected (right) images in the GALEX FUV {top) and NUV (bottom) filters. 



FUV and NUV images are distributed all over the galac- 
tic disk. Since these FUV and NUV filters are especially 
sensitive to young stellar populations (see Fig. [1]) , this in- 
dicates that most of the recent star formation occurred 
in the galaxy disk, and not in its nuclear region. The 
color composites of Fig. [5] and Fig. [T^] show that the 
extinction-corrected FUV and NUV images are clearly 
tracing stellar populations that appear bluer than other 
regions. The images in U, B, and R (Figs. [16] and flT)) 
show the transition of the dominant emission from young 
stellar populations in the SF-regions to older stellar pop- 
ulations in the bulge and center of the galaxy (Fig. [5^) . 



While the extinction corrected image in U is still sensi- 
tive mostly to young stellar populations, the i?-band im- 
age shows that the older stellar populations in the small 
bulge of NGC 959 become the dominant light source at 
redder wavelengths, as expected. 

Among the FUV-i? filters, the largest morphological 
change as a result of the dust-extinction correction occurs 
in the optical images, especially in the i3-band image. 
In the FUV, NUV, and U , the light from the younger 
stellar populations dominates even without applying our 
dust-extinction correction. Correcting for dust extinc- 
tion therefore strengthens already discernable galactic 
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Fig. 16. — Comparison of observed (left) and extinction-corrected (right) images in the VATT U (top) and B (bottom) filters at GALEX 
resolution. 



structures, but does n ot dramatically c hange the appar- 
ent morphology (cf., [Windhorst et al.l [2002). The R- 
band image, on the other hand, is dominated by light 
from older stellar populations and suffers to a lesser de- 
gree from the effects of dust. The S-band samples both 
younger and older stellar populations (see Fig.[T]) and suf- 
fers a larger dust extinction than V (see Fig. [Tj] and Ta- 
ble 2). After correction for extinction, light from younger 
stellar populations that is initially largely hidden behind 
the dust becomes visible, causing the galaxy morphology 
to change relatively more in B than in other filters. 
Previous studies have compared the morpho- 
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bands 

1999; Block et al 



appe arance in B- and N IR 
( e.g.. iBlock fc Pueraril [l99l 



Eskri dge et al.l 120021; 
[200^. [Eskridgc et al 
logical classifications in 
galaxies from the Ohio 
Spiral Galaxy Survey, 
correlation between 
and NIR (see their 
the H-haiid 



120011: 



H- or Ks- 
Elmegreen et al.l 



Buta fc Blo^ 

IBlock et al. 20041; ISeigar" 



2001 



et al.l 

(2002) compared the morpho- 
B and H for ~200 spiral 
State University (OSU) Bright 
They found a relatively good 
the classifications in the optical 
Fig. 2 and Table 2). On average, 
classification was found to be ~ 1 T-type 



earlier than the optical one. Other studies, however, 
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Fig. 17. — Comparison of observed (left) and extinction-corrected (right) images in the VATT R filter at GALEX resolution. 



found n o corr elation between optical Hubble type 
(jHubbld 119261 ) and dust-penetrated morphological 
classes for smaller samples of 14-36 galaxies o bserved 
m K , (e.g.. iB fock & Pucrari 1999; Elmcgrcc n et al 
19991: iBlocket al..,2q01.: ,Buta fc Block ,2001: Bloc k et al 
2004t ISeigar et al.l 120051 1 Since the Hubble types 
are based on blue photographic plates, the Hubble 
classification can be sig nificant ly af fected by dust ex- 
tinction (e.g.. Iwindhorst et al.l 120021 ). i^T^-band images 
suffe r only 10% of the extin ction in V (Fig. [T3] and, 
e.g., iMartin fc Whittetl I199O0 and therefore show the 
stellar distribution with much smaller effects from dust. 
Therefore, if the amount of dust in (and in front of) 
a galaxy is significant, the optical _B-band is affected 
accordingly, possibly resulting in drastically different 
morphology than suggested by the NIR classification. 
As this was seen in only a small subs et of the OSU 
Brigh t Spiral Galaxy Survey samples (jEskridge et al.l 
l2002t ). the galaxy samples u sed in the if^-band studies 
by, e.g., iBlock fc Pueraril ()1999f l may have selected 
dustier galaxies. 

The apparent B-band morphology of NGC 959 be- 
fore and after extinction correction (bottom panels of 
Fig. [16]) does not differ as drastically as the optical 
versus Kg morphologies reported in the above studies. 
But some regions, such as the SE side of the bulge at 
[x, y\ ~ [142, 148], become much more prominent after ex- 
tinction correction and the distribution of light becomes 
similar to that seen at 3.6 and 4.5 /xm (see bottom panels 
of Fig. [T3| . This means that if a galaxy contains a large 
amount of dust along major structures, such as a bar or 
spiral arms, it is possible that the _B-band morphology 
can change significantly after correction for dust. 

5. DISCUSSION 

Correction for dust extinction is an important, yet 
challenging issue when studying stellar populations in 



galaxies, because exti nction has a similar effect as stellar 
population age (eg., Gordo n et al.l 11 9971 and metallic- 
ity (e.g., Worthey 1994; Ka virai et a l. 2007b). Failure 
to correct for extinction will render analyses of stellar 
populations highly uncertain. Commonly, the FIR/UV 
flux ratio is used to measure "the" extinction within 
a galaxy, given as a single averaged extinc tion value 
(e.g. [Buat & Xu 1996; Calzctti ct al. 2000; Bo selh et al] 
120031: iKong et al.l 120041) . or as a onc-dimcnsional (i.e., 
azimuthally averaged) ra dial extinction profile (e.g., 
iBoissier et all 12004 [20Q5I ). This method assumes that 
the distribution of dust is relatively uniform across the 
entire galaxy, or representable by a simple radial ex- 
tinction gradient. In reality, the distribution of dust 
is complex, following galactic structures such as SF- 
regions, bars, and spiral arms. It may also be affected 
by nearby companions or satellite galaxies. Therefore, 
a more detailed analysis of the spatial distribution of 
dust is needed. To obtain a two-dimension al distribu- 
tion of dust extinction, so me studies (e.g., Sc oville et al.l 
120011: iCalzetti et alll2005[ ) use ratios of Hydro gen recom- 
bination lines such as H/3/H7, Ha/H/3, or Ha/P aa, that 
have known intrinsic values (jOsterbrockl 119891 ) . This 
method is applied to measure dust extinction in H li re- 
gions within some of the nearest galaxies, but is appli- 
cable over only a small fraction of an entire galaxy. Our 
method, based on the observed I^-to-3.6/xm ratio and 
model SEDs, is able to map the full two-dimensional 
distribution of dust. Fig. [10] and Fig. [TT] demonstrate 
that the actual distribution of extinction within NGC 959 
does not follow a simple radial trend and that dust is not 
concentrated only in the most actively star-forming re- 
gions; it is present throughout, tracing the complicated 
structures from the galaxy center all the way to the outer 
regions of the disk, where the S /N in the images becomes 
the limiting factor. 
At this point, we would have liked to compare 
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Fig. 18. — a) Radial extinction profiles for the six nearby spiral galaxies from'Boissier ct al.' ("2004), with distances expressed in units of 
^25- While some galaxies show a bump or an upturn in the outermost regions of their disk, the extinction generally decreases from the center 
of a galaxy outward. 6) Visual extinction, Ay, as a function of distance from the center of NGC959. Each dot represents the extinction 
measured in a single pixel. Averages within 0.1 kpc annuli are plotted as asterisks when including pixels with Ay = magarcsec"^, and 
as open diamonds when excluding such pixels. The error bars represent the standard deviation within each annulus. Twelve pixels with 
Ay ^0.5 magarcsec"-^ and >1.3kpc {R}^0.6 R25), are not shown here. Note that the very outskirts of NGC959 may suffer from the 
light-blending from pixels with S/N<3 in one or more of the filters, hence rendering the Ay measurements more uncertain. 



our results for NGC 959 with exti nction m e asurem ents 
from other r nethod s. While Esip ov et alJ ()1991[ ) and 
I Taylor et al.1 (|20G5) have studied this galaxy, they did 
not analyze the internal dust extinction, so we cannot 
directly compare our results to prior work. However, ra- 
dial extinction profiles have been anal yzed in other galax- 
ies (e.g., Q anscn ct al. f 994; Boissi cr et al l [2004. 2005- 
iCalzetti et al.ll2005tlHolwerda et al.ll2009( ). iJausen et all 
(|1994D used the special geometry offered by a nearly 
edge-on disk and a large bulge to demonstrate that the 
maximum extinction in the dust lanes of two galaxies 
(UGC 3214 and UGC 3065) decreased outw a rd wit h dis- 
tance from the minor axis. iBoissier et all (|2004f l used 
the azimuthally averaged FIR/UV flux ratio to measure 
radial extinction profiles for six nearby late-type spiral 
galaxies. Their results are reproduced in Fig. [TSk . where 
radius is expressed in units of R25, the major axis ra- 
dius at the ruB = 25.0 magarcsec"^ isophote as listed in 
the RC3. Each galaxy shows a general trend of decreas- 
ing extincti on from the c enter to the outer regions of 
the galaxy. ^B oissier et al.l ([2005) studied the radial pro- 
files of extinction in M 83 using different methods — the 
Total-IR/FUV luminosity ratio, the UV spectral slope, 
and the Balmer decrement. Their Fig. 2 of shows that all 
three methods give similar results: a general decrease of 
extinction with radi us, with a smal l uptu rn at the outer 
edge of M83's disk. ICalzetti et all 12005) measured the 
Ha/Paa ratio in Hii regions to study the distribution of 
dust extinction in M 51. Their Fig. 14 shows the distribu- 
tion of individual extinction measurements as a function 
of radial distance from the center of M51. While there 
are some relatively highly extincted H 11 regions at outer 



radii — corresponding to the bump at i?/i?25~0.8 in 
the profile for M51 in our Fig. [T8k — the distribution 
docs follow the general decr easing trend wit h increasing 
radius. On the other hand, iHolwerda et al.l (|2009ll used 
an occulting galaxy pair to measure the distribution and 
amount of dust via the optical depth against the back- 
ground galaxy, providing measurements that extend be- 
yond the optically visible disk of the foreground galaxy 
(see their Fig. 2). They showed that large amounts of 
dust can exist even in the outermost parts of spiral galax- 
ies, where these dust features are usually undetectable by 
other means (see their Figs. 11 and 12). 

If we see a similar radial trend for our estimated 
dust extinction in NGC 959 as in the studies above, it 
would lend additional credence to the reliability of our 
method. Fig. [T8b shows the radial distribution of Ay 
for each pixel in NGC 959. The galaxy center is lo- 
cated at [x, 2/] = [151, 148], and the radius is expressed 
in units of R25 ■ Also indicated are mean extinction val- 
ues computed within 0.1 kpc bins in radius, with (aster- 
isks) and without (open diamonds) including pixels with 
Av = magarcsec"^. The error bars represent stan- 
dard deviations for the distribution of Ay values in each 
such bin. While Ay for individual pixels spans a wide 
range at each radius, the azimuthally averaged extinc- 
tion, Av.-mca,n{R), clcarly decreases from the center out- 
ward. At the outermost radial bin, however, the average 
extinction shoots up to Ay,mcan(-R) > 0.3. This is caused 
mostly by the pixels with high Ay within the "interest- 
ing regions" discussed in §3.6.3 (see also Fig. (TU]), which 
may or may not be associated with NGC 959. The larger 
uncertainties for the outermost bins also reflect the fact 
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that these pixels only marginally exceed our minimum 
S/N criterion in one or more of the filters, hence possibly 
affecting the reliability of (3v and Ay- Yet, some pixels 
with relatively high Ay (0.2 < Ay < 0.5) indicate the ex- 
istence of dust in the outermost regions o f NGC 9 5 9, per - 
haps analogous to the findings of Holwer da et al.l (|2009f ) . 
The extensive tests described in previous sections, as well 
as the general agreement with results from other studies 
— although they involved different galaxies — give us 
confidence that our method produces reliable measure- 
ments of the spatial distribution of extinction by dust 
within a galaxy. 

An important lesson from the application of our 
method to NGC 959 is that we are able to map the two- 
dimensional distribution of dust extinction even from the 
low-resolution images, which show no conspicuous dust 
features silhouetted against the stars in Fig. [5b. Our 
initial expectation was that the extinction map created 
with our method might be relatively featureless, with 
specific dust features smoothed out. Instead, a map with 
a complicated pattern that follows the galactic structure 
emerged (see Fig. [TU|) . A similar, but more detailed dust 
distribution is recovered when the analysis is repeated on 
images at the higher IRAC 3.6 /im resolution (Fig. [TT|) . 
Pixels with large Ay in Figs. [TUl and [TT] clearly trace the 
location of silhouetted dust features in the higher reso- 
lution VATT (Fig. Ell) or HST (Fig. [H]) color compos- 
ites of NGC 959. These results stress that the contrast 
between regions with high and relatively low extinction 
can be large within a galaxy disk. At resolutions of 2'.'2 
{Spitzer /IRAC) and 5'.'3 {GALEX) or ~110-250pc at 
the distance of NGC 959, our method is able to reliably 
and meaningfully generate a two-dimensional distribu- 
tion of the dust extinction. Such modest resolutions are 
accessible and our method will be applicable in galaxies 
well beyond the Local Group, where the individual dust 
features may no longer necessarily be spatially resolved. 

Our method will also be useful to study dust extinction 
in galaxies at much larger distances. Once HST /WFC3 
and JWST are in operation, galaxies at z~ 0.3-0.4 will 
be easily accessible for such studies, and galaxies at even 
larger redshifts might be reachable. WFC3^^ has two 
imaging channels: UVIS covers 0.2-1.0 /itm with a pixel- 
scale of 0'.'04 pixel"^, while NIR covers 0.9-1.7/im at 
0'.'13 pixel"^. JWST^"^ will also have two imaging cam- 
eras: NIRCam, covering 0.6-5 /xm at 0'.'032 or 0'.'065 
pixel"^ (Short versus Long Wavelength Channel), and 
MIRI, which covers 5-27^m at ~0'.'ll pixel'^ At 
z>0.3, the rest-frame V and 3.6 /im bands shift pro- 
gressively further into the near- to mid-IR regime cov- 
ered by HST /WFC3 and JWST, allowing one to apply 
our method to large numbers of intermediate redshift 
galaxies. The only limiting factor is the apparent size of 
the galaxy at these redshi fts. Up to z ~ 1.6, the angu- 
lar scale becomes smaller (jWrightl [20061 ) and each pixel 
samples a larger surface area within a galaxy. At z ~ 0.4, 
the central wavelengths of the V and 3.6 /im bands shift 
to 0.77 ^m and 5.04 ^m, where UVIS and MIRI are ex- 

Space Telescope Science Institute (STScI), Wide Field Cam- 
era 3: 

http : //www . stsci . edu/hst/wf c3 

1"' NASA, The James Web Space Telescope: 
http : //www. jwst .nasa. gov 



pected to deliver resolutions of O'.'OS and 0'.'195 FWHM, 
respectively. This corresponds to '^420pc and '^Ikpc, 
where the latter sets the relevant resolution for our 
method. While dust lanes will certainly not be visible, 
we expect our method to still produce meaningful maps 
of the variations in extinction on scales of ^ 1 kpc, as 
long as a galaxy is at least several kpc in diameter. 

Another a pplication of our rn ethod is to investigate 
the result of iTavlor et al.l (|2005f) . that the color in the 
outer regions becomes redder in the majority of late- 
type spiral and irregular galaxies. This reddening may 
be caused either by a change in stellar populations or 
be due to the presence of dust, which is usually only 
apparent in higher resolution images when silhouetted 
against a relatively b right stellar background (as in, e.g., 
iHolwerda et al.ll2009[ ). In NGC 959, a nearly face-on late- 
type spiral galaxy, we found evidence for the existence of 
a moderate amount of extinction in the outermost regions 
of its disk (Fig. fTU]) . Our method will also be useful to 
study the dust content of elliptical and lenticular galax- 
ies. In recent years, large amounts of dust were discov- 
ered to exist in elliptical and lenticular galaxies, as well 
as in the halo of spiral galaxies (c.g^, K ayirai et al. 2007a|; 
lOosterioo et al1l2007al fg: lEmonts et aLil2008i r Since very 
little SF-activity is ongoing in these galaxies, methods 
based on the FIR/UV flux ratio or on Hydrogen recom- 
bination lines are not as useful as they are for actively 
star-forming galaxies. Our method, which only depends 
on images in V and 3.6 /im filters, is well-suited for a 
study of the distribution of dust within early-type galax- 
ies. We will present a more detailed study of dust distri- 
bution in spiral galaxies and in a small number of ellip- 
tical and lenticular galaxies — a total of 45 galaxies of 
all types — in a subsequent paper (Tamura et al. 2009c, 
in preparation). 

A potential future use of our method is to estimate 
the distribution of dust in simulated galaxy models (e.g., 
iCroton et al.ll2006l ). A detailed two-dimensional analysis 
of the dust extinction in a large number of galaxies would 
help understand the properties of dust for different types 
of galaxies or galactic structures, such as SF-regions, 
bars, spiral arms, and inter-arm regions. Current sim- 
ulations are able to model the extinction-free SEDs for 
many galaxies, while the treatment of internal dust ex- 
tinction is still broad-stroke. For such models, it will 
be useful to construct a database of galaxies of different 
morphological type and mass, to constrain age, metallic- 
ity, and the amount and spatial distribution of the dust. 

6. CONCLUSIONS 

In this paper, we presented a new method for estimat- 
ing the extinction by dust within galaxies by comparing 
the observed V-to-3.6 iim fiux ratio, (3y, to theoretical 
SED models. Using a pixel-based analysis, our method is 
able to estimate the two-dimensional distribution of dust 
extinction within a galaxy. As a proof of concept, we ap- 
plied this method to NGC 959, a nearby late- type spiral 
galaxy. From a pCMD, constructed using an additional 
C/-band image, we robustly selected pixels dominated ei- 
ther by the light from younger stellar populations, or 
from older ones. Since their intrinsic y-to-3.6 /zm flux 
ratios differ, they were treated separately in our anal- 
ysis. We presented a two-dimensional map of the vi- 
sual extinction. Ay, that closely resembles the observed 
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distribution of SF-regions and underlying galactic struc- 
tures (including a newly identified bar) , and which traces 
the distribution of 8.0 /xm PAH emission. Although dust 
features are inevitably smoothed out to some extent due 
to light-blending and the low spatial resolution of the 
images, we were able to construct a two-dimensional ex- 
tinction map with sufficient detail to delineate the struc- 
ture of dust features within the disk of NGC 959. We 
then presented original and extinction-corrected views of 
NGC 959 from the FUV through MIR. Through a series 
of tests, we demonstrated the validity of our results and 
method. 

Our method has several advantages over other meth- 
ods based on, e.g., the FIR/UV flux ratio, UV spec- 
tral slope, or Hydrogen recombination line ratios. At 
its core, our method only depends on images in two 
relatively common broadband filters, V and 3.6 /xm (L- 
band) , and is therefore applicable continuously across the 
face of a galaxy. We exploit the fact that the wavelength- 
dependent extinction by interstellar dust reaches a mini- 
mum near 3.6 /im while it increases toward shorter wave- 
lengths. In y-band, we are sensitive to dust extinction, 
but fairly insensitive to age and metallicity effects com- 
pared to UV-blue filters. We demonstrated that the in- 
trinsic V-to-3.6 iim flux ratio, f3v.o, is well-behaved over 
a wide range in stellar age and metallicity. While fJv.o 
depends on age more strongly than on metallicity, this ra- 
tio stays relatively constant for older {t > 500 Myr) stellar 
populations, and occupies a relatively narrow range for 
younger stellar populations. As a result, we can simply 
compare the observed and intrinsic y-to-3.6/im flux ra- 
tios (after taking mixing or superposition of stellar pop- 
ulations into account) to estimate the amount of dust 
extinction. Ay, in each pixel. To translate Ay to bluer 
filters, knowledge of the metallicity become more impor- 
tant (or, alternatively, the uncertainty increases) due to 



the metallicity-dependence of the extinction curve. This 
simplicity allows our method which is mostly auto- 
mated with only a few manual parameter-settings for 
each galaxy — to be applied to a large number of galaxies 
in a very short time. 

Since our method does not require visual confirmation 
or identification of individual dust features, it is appli- 
cable to any galaxy beyond the Local Groups, if rest- 
frame V and 3.6 /xm images are available with at least 
several resolution elements across that galaxy. This of- 
fers the possibility of applying our method to HST/ ACS, 
F5T/WFC3, and JWST NIRCam and MIRI images to 
study the two-dimensional distribution of dust not only 
in the local universe, but also for higher redshift galaxies. 
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